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At high e-foldings of expansion, the inflation with the quartic potential exhibits the parametric 
attractor governed by the slowly running Hubble rate. This quasiattractor simplifies the analysis of 
predictions for the inhomogeneity generated by the quantum fluctuations of inflaton. The method 
reveals the connection of inflation e-folding with general parameters of preheating regime in various 
scenarios and observational data. 
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I. INTRODUCTION 



The homogeneous component of inflationary dynamics 
in the simplest case of sing le field fH, 0, 1 3 (see review 
in 0) is usually described in terms of slow-roll approx- 
imation in the evolution of scalar inflaton 0, when one 
neglects both its kinetic energy ^(jP and acceleration (f> 
in the field equations 
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2 ^^r + v{cj)) 



dV 



(1) 



derived from the Hilbert-Einstein action of Friedmann- 
Robertson-Walker metric 



ds^ ^dt^ -a^{t) dr^ 



(2) 



with the scale factor of expansion a{t) ordinary defining 
the Hubble rate H = a/ava terms of its derivative with 
respect to time t as denoted by over-dot 6 = d o /dt. In 
([1]) the Planck mass is usually introduced through the 
Newton constant G by the relation rrif-^ — 1/G, while V 
is the infiaton potential. The slow-rolling implies 
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For the consistency of such the approximation we have 
to estimate the ratios 
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and 
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e < 1, 
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Then, introducing parameter rj by 
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and using we get 



(7) 



Therefore, the slow-roll approximation is substantiated 
at 



e<l, hi < 1. 



(8) 



Constraints ^ are usually satisfied for positive power 
potentials, for instance, at appropriately chosen parame- 
ters and fields (j) > ?npi, which guarantee the infiation of 
Universe with slowly changing Hubble rate with a huge 
total e-folding of scale factor A^tot. = InOcnd/ainit., where 
subscripts stand for the moments of ending and initiating 
the inflation, correspondingly. 

Note, that the above speculations have not involved 
any initial data on the field evolution. This can be mean- 
ingful, if only the evolution possesses properties of attrac- 
tors, which allow a stable dynamical behavior, that soon 
forgets about a start. 

The quasiattractor for the infiationary dynamics with 



the quadratic potential V ■ 



was found in ^ , which 



modernized the consideration of dependence on the ini- 
tial conditions in the homogeneous case as was done in 
pioneering papers of 0, So, the phase space vari- 
ables of dynamical system rapidly tend to a stable point, 
which position depends on the parameter defined by H , 
while the parameter itself slowly evolves during the most 
amount of e-folding in the infiation and it begins sig- 
nificantly to change to the end of infiation at a short 
increment of A^tot., only. Thus, at high total e-folding of 
infiation governed by an appropriate large starting value 
of Hubble rate, the infiation due to the quadratic po- 
tential is well described by the quasiattractor behavior. 
Similar result was obtained for the quartic potential in 
for the homogeneous case, and further attempts to gen- 
eralize the analysis to damping of initial inhomogeneities 
were formulated. 

The analysis of dynamical system for the homogeneous 
infiaton and baryotropic matter was done in ^3] , wherein 
some general features of evolution were derived for an 
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arbitrary potential of inflaton, and several explicit exam- 
ples were studied. 

Another aspect of considering the initial conditions 
concerns for inhomogeneities perturbing both the infla- 
ton and metric. The analysis of inhomogeneity in the 
initial values of scalar field beyond the perturbation the- 
ory as well as relevant references can be found in TT], 
wherein the authors formulated conditions constraining 
the parameters of initial state, when the inflation takes 
place. 

In the present paper we generalize the approach of Q 
to the case of quartic potential V ~ that modernizes 
the consideration in The main goal of our study is 
to show the effectiveness and elegance of quasiattractor 
method to the analysis of inflation dynamics and its con- 
nection with the preheating regime in the case of quartic 
potential in order to draw the conclusion on the con- 
sistency with the observational data versus general con- 
ditions of preheating. The dynamical system relevant 
to the evolution of homogeneous Universe is considered 
in Section |IT1 We find the stable quasiattractor for the 
quartic potential, too. The inflation parameters are in- 
vestigated in Section IIIII wherein they are confronted 
with both the experimental data and predictions of slow- 
roll approximation. Our results are summarized and dis- 
cussed in Conclusion. 



II. DYNAMICAL SYSTEM 



The action of inflaton is of the form 



with the quartic potential 



(9) 



(10) 



In the Friedmann-Robertson-Walker metric it leads to 
the evolution 



4> ^ -3i?0- A(/)^ 



(11) 



The analysis of dynamical system is simplified, if we in- 
troduce the driving parameter 



/3A 



(15) 



Then, denoting the derivative with respect to the amount 
of e-folding by prime o' = do /dN at = In a — Inoinit., 
we get 

x' = -3x- 2y^z, 
3 



(16) 



-X y + xz, 

3 2 
-X z. 



It is an easy task to confirm, that the system of pB]) 
conserves the Friedmann constraint in the form of (|14p . 

It is important to stress, that z is quasi-constant, if x 
is fixed at a critical point Xc ^ 1, hence. 



•^■init. ^ 



\xlN 



^init. I l + -a;^Ar + 
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so that the monotonic growth of z is slow in comparison 
with the linear increase of e-folding A^ due to a small 
slope proportional to x\ 0, indeed, and one can put 
z « Zinit. until '^x\N ^ 1, i.e. at rather large intervals 
of A^. We will see later, in fact, that the stable critical 
point x\ scales as x\ « /3/(A^tot. — A^) at A^tot. — A^ > 1 



and /3 



so that the growth takes the form 
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under the Friedmann equation 



Therefore, at A^ ^ A^tot. the increase is quite linear, but 
the coefficient of Zinit. is actually suppressed Zinit. <C 1 
by the conditions of attractor, as we will consistently see 
later. 

Thus, we can treat the phase-space evolution in the 
{x, plane as the autonomous dynamical system with 
external parameter z, when the driftage of z gives rise to 
the sub-leading correction. 
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(12) 



at = SttG. 

Following the method of , let us introduce the phase 
space variables 



4/ A K0 



related due to (|12p by the constraint 



(13) 



x^+y^^l. (14) 



A. Critical points 

Solving a;' = = at X 7^ and y 7^ 0, we get the 
relation^ 



ixy — — 2z, 



(19) 



^ The critical point x = 
satisfy constraint I I14I I. 



is off interest, since it does not 
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while 



B. Stability analysis 



/(x-^) = - + l-z] =0 



(20) 



The cubic polynomial f{u) in (j20p produces two extremal 
points, since its derivative is equal to 



du 



u{3u — 2), 



hence, Ui — 0, U2 ~ ^- Therefore, / has got positive 
roots, if only 



f{u2) < 0, 



that yields 



(21) 



Constraint ((2T|) guarantees the existence of critical points 
Substituting 



z^^-sin^x, at < X ^ ^, 



(22) 



allows us to write the real solutions of ((20|) in the follow- 
ing form 



1 



2 2 



a:,^--^l + 2cos( -x + -7r^) !>, ^ = {0,±1}. (23) 



At £ — +1 quantity ^ 0, which is not a physical 
point, at all, while at £ = — 1 and £ = 



(24) 



correspondingly. Therefore, the case of ^ = is covered 
by the case oi £ — —1, if we expand the interval of x to 

< X < TT, (25) 

that allows us to put the critical point^ to Xc = x^. 
Due to (fT9|) . we can write down 

fix^) = xHx' + y'^l), 

so that the critical points satisfy the Friedmann con- 
straint of (full. 



^ One could change the convention by setting — tt < x < and 
Xc = XQ, but we prefer for l|25|l . when the case of interest with 
<C 1 is spectacularly reached at x ^ 0. 




At and y — Uc + y, the evolution linear in 

{x, y} reads off as 



(26) 



wherein we have used (fT9|) . 

The linearized Friedmann constraint in (I14|) gives 

= -22/3y, (27) 
so that is reduced to the single equation 



y 



9/2 2 

xt 



Therefore, at 



xl < 



3' 



(28) 



(29) 



the deviations {x, y} will be damped, while at Xc — > 0, 
the fluctuations will decline as 

The same fact can be found in a general manner. So, 
the eigenvalues of matrix in are the following: 



A, 



6xi 



v+ 



while the corresponding eigenvectors are given by 
However, the condition of (I27p can be rewritten as 



y 



(30) 



(31) 



(32) 



Therefore, the eigenvector v+ is irrelevant to the consid- 
eration, since it drives out of constraint (|27p . while V- 
is consistent with (|27p . and the critical point is stable at 
A- < 0, i.e. under (^5)) . 

Thus, we have established the quasiattractor for the 
quartic potential of inflaton. The stability of critical 
point is controlled by smallness of z, that means large 
enough value of initial Hubble rate. 



III. INFLATION PARAMETERS 

Staring the evolution of Universe filled by the inflaton 
from a high initial Hubble rate, i.e. at Zjnit. 0, but 
i/j^ij < TTip]^, causes a rapid entering the quasiattractor 
regime, when the expansion is described by the slow drift 
of Hubble rate. It means that a small domain of space 
is inflationary increased to a huge size, and the inflaton 
can be approximated by an almost homogeneous field 
versus the spatial coordinates, while its value and time 
derivative evolve in accordance with the drift of critical 
points. 
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A. Homogeneous limit 



The quasiattractor is the reason for the cosmological 
evolution forgets its primary origin to the leading homo- 
geneous approximation. Nevertheless, one could consider 
some properties of attractor dynamics close to the end of 
inflation, since it further enters the properties of observed 
inhomogeneity. 



1. Acceleration conditions 

The accelerated expansion takes place at a > 0, hence, 
a 



a 

that is reduced to 



^2 2ij2 



3x^ > -1, 



(33) 



wherein we have used Therefore, the Universe fol- 

lows the inflation at 



2 , 2 
"^c '''c.cnd 



(34) 



i.e. in the region, where the quasiattractor is in action, 
indeed. 

The accelerated regime ends at 



2 
3' 



:,cnd ~ \ 2 ^c,cnd2/c,cnd 
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(35) 



equivalent to Xcnd = f and 



48_A 
1? 



The inflation condition of (|34p is supplemented by ac- 
companiments 



(37) 



Note, that to the end of inflation the field takes the 
value of the order of Planck mass independently of cou- 
pling constant A, while the corresponding squares of Hub- 
ble rate and kinetic energy are scaled linearly in the cou- 
pling constant, and they are actually suppressed to the 
appropriate Planck mass powers at A ^ 1. 



2. Amount of e-folding 

Once the attractor is in action, the total amount of 
e-folding during the inflation can be obtained by inte- 
gration of z-component in the system of (jl6p at x = Xc, 



Zend 

2 r dz 



(38) 



In the limit of z-mit. ^ 0, we deduce (1381) by transforma- 
tion to the integration versus x, so that 



7r/4 



cot X dx 



Xinit 



approximated by 



(0) 



f -f 2cos(|x- f tt) 
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\/3 /■ dx 



(39) 



2 J X- 

Xinit. 



(40) 



to the leading order at Xinit ~* 0- T^^^ consequent ex- 
pansion in (|39p gives 



(0) 



(1) _ 
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TT 
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^ '^init. 





1, 4x,„,, 
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(41) 



The inversion of (|4ip results, for instance, in 



X 



(0) _ 



2t:N, 



(0) 
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(42) 



To the leading order at A'tot. S> 1 we deduce the initial 
data 



1 



1 1 
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(43) 



while j/c.init. = 1 ^ ^c.init.- Notc, that the initial data 
for the field and its velocity {y, x} mean their values just 
after entering the quasiattractor^ while the actual data 
are irrelevant, since the quasiattractor rapidly adjust the 
field and its velocity to in accordance to the initial 
Hubble rate, i.e. Zinit.- 



B. Inhomogeneity 

Quantum fluctuations of infiaton near its homogeneous 
classical value result in the primordial spatial perturba- 
tions of energy density, involving the scalar and tensor 
components of spectrum versus the wave vector k at the 
moment, when the fluctuation comes back from the hori- 
zon, i.e. at fc = a{t) H as following: 



(44) 
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The spectra can be accurately evaluated in terms of 
quasiattractor dynamics by 



87r2 z^xl ' 



-r, /,^ 6A 1 



(45) 



where z corresponds to the Hubble rate H, and it is ex- 
pressed by amount of e-folding N between the moments 
of horizon exit and inflation end, as that can be found in 
the right analog of in the leading order at TV ^ 1. 
So, we redefine N = Inocnd — In a. Then, 



(46) 



which are consistent with the slow-roll approximation for 
the inflation with the quartic potential. The sub-leading 
corrections are given by the substitution following from 



1 

TV 



27r 



2ttN + A^/i' 



The ratio 



Vt 9 

r = :^=48x2«l, 



(47) 



(48) 



determines the relative contribution of tensor spectrum. 
In the leading order 



16 

iV' 



(49) 



in accordance with the slow-roll approximation, again. 
The spectral index is defined by 



ns - 1 



dlnT's 



d In A; 

It can be calculated under the condition^ 



(50) 



ln-^ = -iV-21n 



, (51) 



so that 



ns - 1 



(l-3a;2)(2-3x2)' 



(52) 



that is reduced to the slow-roll result in the leading ap- 
proximation 



ns - 1 



^ N 



(53) 



^ Note, that e-folding counts backward, from the end of inflation 
to the moment t, hence, the derivatives in I I16I I change sign. 



To the same approximation 

na — 1 ~ r. 

^ 16 

The running of spectral index is given by 

dns ^{6-16xl + 9xi){l- xl) 



d In A; 



-36x^ 



(2-3x?)3(l-3x2)3 



approximately equal to 



-27xt 
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(54) 



(55) 



(56) 



As it was noted in the evaluation of a in the 
framework of slow-roll approximation is quite compli- 
cated, while the quasiattractor dynamics allows us to get 
it straightforwardly by means of taking the derivative. 
Though one could use (|?T|) to the leading order, when 
dk/dN ~ — 1, that straightforwardly yields the spectral 
index as well as its slope derived above. 

Thus, the quasiattractor behavior provides us with the 
direct way of calculating the properties of primary in- 
homogeneity in terms of two parameters: the coupling 
constant A and number of e-folding N relative to the ob- 
servational scale k. 



1. Constraining a maximal N 

Following the method elaborated in [l^ , let us restrict 
the scale k = aH relative to the observation of inhomo- 
geneity, as it comes from the both known and suggested 
history of Universe evolution by considering the ratio 

k aH 



qoHq aoHo 

where Hq is the present-date Hubble rate 

i7o = 1-75 • 10~^i/iTOpi, /i~0.7, (57) 

while flQ is the present-date scale factor, that can be set 
to flo = 1 by definition of length unit. We decompose the 
evolution in several consequent stages: the inflation, re- 
heating, radiation dominance up to the epoch of equality 
with nonrelativistic matter, matter dominance, that give 

^ ^ ^cnd ^rch. ^cq. 

flO lend arch. Ooq. ^0 



H _ H Hcnd Hcq. 



(58) 



^^0 Hcnd Hen. Ho 



The data yield 



Hn 



5.25-lO'^h^flliHo, "''q -|^°q- ^ 219 /t»M, (59) 



where VIm is the fraction of matter energy. The quasiat- 
tractor gives 



-"end 



'end 



(60) 
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The scale ratios are related with the ratios of Hubble 
rates in two subsequent stages "6" "c" by means of 



7j2 tt2 



3(l+iuc 



(61) 



where Wc is the state parameter at the stage "c" , deter- 
mined by the ratio of pressure to the energy density. For 
the radiation stage we set Wcq. — |. 
Then, we arrive to 



= - In 



— |- +ln219f^M/i+ Jln^ 

dQllQ Z O 



IniV 



1 



3(1 + Wroh.) 
1 

3(1 + Wroh.) 



tt2 

In^ 



In ^-'^ 



(62) 



1 , H 

7^^^ — 
4 mpi 



Let us parameterize the Hubble rate at the end of reheat- 
ing by a scale /irch., so that 



Hi 



3mp]^ 



Mrch. 



while 



Therefore, 



tt2 
^end 

ml. 
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Vs. 



(63) 



(64) 



iV = - In ■ 
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ln219f7M/i+ Trln- 



1 



1 + Wroh. 
1 



In AT 



3(l+Wroh.)/ 2 



(65) 



1 



3{l + Wrch.)J 3 m 



In ^ ^'■'=11 ■ 



Note, that (|65p is independent of current matter density 

Astronomical observations are naturally restricted by 
the distances less than the event horizon given by the 
inverse Hubble rate of present day Hq. Therefore, the 
maximal value of e-folding A^o is reached at the wave 
number k kg = uqEq, when (j65p contains two pa- 
rameters determined by the mechanism of reheating. In 
addition, experimental data give the spectral density (l3 | 



2.5-10" 



(66) 



at the scale k ~ 0.05 — 0.002 Mpc ^ greater than fcg, 
nevertheless, this value of Vs rather slowly evolves with 
k. 

The dependence of Nq versus the scale of reheating 
firch. and state parameter Wrch. is shown in Fig. [1] for 



generic cases of stiff matter w = 1, radiation w = and 
dust w = 0. It is spectacular, that, first, the dependence 
versus the scale is well approximated by logarithmic law, 
second, all these reheating regimes converge to a single 
point at /2rch. — 0.33 • 10^^ GeV, and third, the radiation 
regime gives the maximal amount of e-folding'* equal to 
Ao.rad. ~ 64 independent of the scale characterizing the 
reheating stage. The reheating scale /lioh. for the inter- 
section of curves with different values of state parameter 
Wroh. in Fig- fflis given by the condition of nullifying the 
coefficient in front of factor 1/(1-1- Wroh.) in (pS]). so that 



Prch. 



9 Vs 



16 A^o'rad. 



(67) 



Our numerical calculations exhibit the oscillations of 
the field around the minimum of quartic potential with 
Wreh. ~ 1. However, this value can be modified by the 
both corrections to the potential near the minimum and 
inflaton interaction with ordinary fields. 



90 



lO" 10"' 

fiich., GeV 



FIG. 1: The maximal number of e-folding A'o versus the re- 
heating scale in different regimes of reheating state parameter 



Wroh. = 1 (long-dashed line), Wr, 



I (solid line), Wrah. ~ 



(dashed line), w^roh. 
(short-dashed line). 



oo (dotted line), and Wreh. 



We have also shown two marginal cases. The first case 
corresponds to Wroh. = — 5 and it could appear, if we sug- 
gest the inflaton oscillations near the potential minimum 
neglecting the gravitational damping during the period. 
Then, the virial theorem gives averages of kinetic and 
potential energies as (a;^) = ^(y^), that reproduce the 
mentioned state parameter. 

The second case corresponds to a quite instantaneous 
falling of Hubble rate to its reheating scale, so that 
Wich. 00. Such the regime could happen due to a 
specific tachyonic preheating fT3| with a rapid decay of 
classical inflaton to quanta, for instance. However, the 



The value of A^o.rad. is obtained by numerical solving l|65| l at 
u. = i. 
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jump of Hubble rate is the indication of energy jump, 
that is unreahstic, at all. Nevertheless, the situation with 
w > 1 takes place, if the potential minimum is shifted to 
a negative value KninO < 0, that gives 

1 ^2 _ ^ 1 02 _ y^^,^ 

W = l-> > 1. 

y^+V \4>-^+ Knin 

Therefore, the kinetic energy and Hubble rate will rapidly 
fall down during short increment of e-folding, while the 
state parameter will follow 3> 1, until the potential 
will become to grow. Then, we suggest the existence 
of second local minimum of V at V^(O) = giving the 
flat vacuum. The decay of flat vacuum to the Anti-de 
Sitter point Knin < is preserved by the gravitational 
effects [H, [l^, so that the flat vacuum can be stable. 
The barrier between two minima of potential suggests the 
presence of negative second derivative of potential with 
respect to the field, that can switch on the preheating 
mechanism at the scale of potential barrier denoted by 
, which could be much less than the energy density at 
the end of inflation. Thus, in Fig. [T]thc region between 
the dotted and long-dashed lines can be actual in the 
extended version of quartic inflation at |Knin| <C ra^-^. 
Though, we cannot point to any verified realistic model 
of such the scenario to the moment. 

The observational data deal with the wave number k 
about two orders of magnitude less than /cqj that diminish 
the corresponding number of c- folding by the value about 
(5iV ~ 4 - 5. 



2. Comparing with data 



X ■ 10" 
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FIG. 2: The experimentally admissible region in the plain of 
{N,X) with la and 2a contours as follows from (|68p . (|69|l . 
The shaded region is allowed after taking into account the 
correlation of spectral index ns with the fraction of tensor 
perturbations in the energy density r at 2a level. 

fraction of tensor perturbations in the energy density r, 
that is shown in Fig. [3l The correlations restricts the re- 
gions of ns admissible for the quartic-potential inflation, 
as shown in Fig. [2] by the shaded domain. 



For the flat Universe with cosmological constant, re- 
cent WMAP data [ll] at fc = 0.002 Mpc^^ gave quite 
precise values of 

ns = 0.95lt°:°i^, Vs - 2Mt"o fe ' 10"', (68) 

while the slope of ng is not so restrictive 

a = -0.102i^:Slg. (69) 

Then, we extract the amount of appropriate e-folding 

N = 6ltll (70) 

that yields the coupling constant equal to 

A= 1.6t?:^ 10"^^ (71) 

which follows from Fig. [2] 

Note, that the inflation predicts the significant sup- 
pression of slope a with respect to the central value in 
(|69p . Therefore, an improvement of accuracy for the mea- 
surement of slope could serve for discriminating the in- 
flation models. 

Further, the WMAP data gave a strong correlations 
between the spectral index of scalar perturbations ng and 



r 




0.9 0.95 1 1.05 



FIG. 3: The WMAP data on the correlations between the 
spectral index of scalar perturbations ns and fraction of tensor 
perturbations in the energy density r, as represented by la 
and 2a contours. The line gives the predictions of quartic 
inflation versus large amount of e-folding A'^. The dot stands 
a,t N = 79, while the rectangle does at A = 59. 

Thus, the A0'*-inflations could be marginally consistent 
with observations, if one suggest the extended version of 
preheating at low scales about /Xj-eh. ~ 10^ GeV with 
the preliminary passing the region with negative values 



8 



of potential as we can conclude from Figs. [T]-[3l An 
example of model for such the potential will be considered 
elsewhere [13] ■ However, to the moment there is no a safe 
realistic model for the preheating scenario consistent with 
the observational data in the quartic inflation. 

IV. CONCLUSION 

We have shown that the inflationary dynamics with 
the quartic potential obeys the parametric quasiattrac- 
tor governed by the Hubble rate slowly evolving with e- 
folding of expansion. The condition of attractor stability 
is preserved by the condition of accelerated expansion. 

The quasiattractor allows us to express the inflationary 
parameters in terms of coupling constant A and amount 
of e-folding N in consistence with the slow-roll approxi- 
mation. Sub-leading terms to the approximation are also 
on hands. 

For the case of quartic potential, we have re-analyzed 
the possible maximal amount of c-folding A^o correspond- 
ing to the scale of astronomical observations measuring 
the inhomogeneities generated by the quantum pertur- 
bations of inflaton just before the end of inflation. It is 



spectacular that at the reheating scale /ircii. 0.3 • 10^^ 
GeV, the value of Nq ~ 64 is independent of the partic- 
ular mechanism of reheating parameterized by the state 
parameter Wrch.- At the low-scale reheating the maximal 
Nq could be increased due to the essential modification 
of potential near the origin. For instance, the field could 
pass the region of negative potential with the further re- 
laxation in the flat minimum at = after overcoming 
the barrier producing the tachyonic preheating. 

Then, the standard quartic inflation with realistic pa- 
rameters of preheating regime is inconsistent with the 
observations of matter-density fluctuations, its spectral 
index of scalar perturbations and fraction of tensor fluc- 
tuations. However, there is the marginal case at the cou- 
pling constant A ~ 6 • 10^^'* and amount of e-folding 
A^ ^ 80 in the modifled preheating at low-scales with the 
negative valley of potential described above, so that the 
quartic potential of inflaton at high flelds with the appro- 
priate modiflcation near the origin is generally still not 
excluded, but a realistic model with such the scenario is 
not known. 

This work is partially supported by the Russian Foun- 
dation for Basic Research, grant 04-02-17530. 
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